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We present a brief survey of the development of nuclear physics towards relativistic 
quark physics. This is followed by a thorough discussion of the quest for the 
observation of the dissolution of nuclear matter into the deconfined quark matter 
(QGP) in relativistic nuclear collisions. Use of strange particle signatures in search 
for QGP is emphasized. 



1 From Nuclear to Quark Physics 

What (we think) we know and what we know we do not know 

The atomic nucleus is a quantum bound state of nucleons (protons and neu- 
trons) each comprising three confined i.e. permanently bound valance quarks. 
The dynamics of nearly massless quarks q — u, d inside a nucleon is fully rela- 
tivistic. Thus the mass of the nucleon and hence of the nucleus, and therefore 
of about 98% of all matter known to us, is understood to be the relativistic 
quark quantum-zero-point energy. Confinement localizes quarks to within a 
volume with radius of about 1 fm in size, about 20-40 times less than (light 
q = u, d) quark Compton wavelength X(q) = hc/m q , with 3 < m q < 15MeV. 
Most important point here is that the indirectly arising interactions, due to 
the properties of the confining strongly interacting vacuum state, determine 
the scale of hadronic size, and thus mass, and thus all the properties of the 
Universe we see around us. Moreover, as we shall argue at the end, there is 
all reason to believe that the hadronic energy scale is a direct descendent the 
unification scale of gauge interactions. 

Research of past 50 years, since the discovery of the pionELcstablished the 
following paradigm of the properties of strong interactions Bu: SU(3) C -QCD 
is the fundamental gauge theory of strong interactions, with quarks q and 
gluons G being the color charged fundamental degrees of freedom; all strongly 
interacting and in our vacuum mobile particles are color neutral. To explain 
the confinement of colored states to localized region of space time the principal 
postulate isQ that the true vacuum state \ V) abhors the color charge; but there 
is an excited state \P) referred to as perturbative vacuum in which colored q, G 
particles are mobile and have their naive perturbative physical relevance. 

The physical difference between \P) and \V) is akin to the difference be- 
tween vapor and ice; in regions of space-time in which 'elementary' vapor 
exists, where \P) replaces |V), a finite energy density the latent heat B has 
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Table 1: Strong interaction vacuum states 
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been deposited. Pushing the vapor and liquid phase similarity we see that 
there may be an intermediate state of the vacuum which we shall call \F). In 
this 'fluid flavor state' the color quantum numbers of quark and gluons are 
mobile, but the dissolution of \V) is not complete, and relatively little latent 
heat is consumed in the transformation \V) — > \F). The principal degrees of 
freedom in the |-F)-vacuum could be effective 'constituent' quarks. 

Thus to dissolve fully the quantum vacuum structure a considerable latent 
(condensation) heat per unit volume needs to be delivered. It is thought that it 
is of the magnitude ,6=0.1-1 GeVfm -3 . The lower value Bp = 0.1 GeVfm~ 3 is 
the energy scale required for the understanding of hadron masses (see Section 
||), the upper value Bp = 1 GeVfm~ 3 arises aside from lattice gauge studies 
of QCDtJ from the necessity to make the deconfining, perturbative state \P) 
inaccessible in low energy reactions, in which free quarks and gluons do not 
become apparent as the fundamental degrees of freedom. Bp is the dimensioned 
quantity which determines the hadronic mass and size scales. The value of 
Bp = 0.1 GeVfm -3 has to be compared to 'normal' matter condensation 
energies which have a magnitude of say O.leV A -3 , a factor 10 24 smaller! 
In order to study states with extreme properties rivaling those found in early 
Universe some 40 [is after its birth we have developed experimental facilities 
in which heaviest nuclei are made to collide at high energies. 

Unclear horizons 

There are many open questions and unresolved puzzles in relativistic quark 
physics. Some are fundamental and border on the comprehension of the origin 
of the standard model, other are simply consequences of the complexity of 
the gauge theory that governs quantum-chromodynamics (QCD), the theory 
of strong interactions. 
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• The fundamental questions are of the type: why do quarks have so many 
different quantum numbers, e.g. i) fractional electrical charge, ii) color 
non-Abelian charge, iii) several (three) families of particles, where (u,d), 
(c, s) and (t, b) are the three doublets of quarks, iv) intrinsic spin. How- 
ever, an elementary object is expected to posses just a few simple prop- 
erties. Are perhaps even quarks not elementary? 

• Understanding of quark confinement riddle falls under the practical cate- 
gory of solving QCD. It is not merely the problem of the vacuum proper- 
ties, but significantly, a problem of hadronic spectrum. The key problem 
that emerges is that we cannot find experimentally many of the bound 
states of quarks that intuitive extrapolations of the current understand- 
ing imply. 

Diverse simple quark bound states other than known mesons and baryons corre- 
sponding to angular and radial quark excitations inside the confining boundary 
are absent. Other expected states, such as the nucleon-anti-nucleon molecule 
(baryonium) which in the quark language consists of the spatially separate two 
quark - two antiquark clusters, qq-qq, have also not been seen. 

The nucleus is made of 3N quarks, but quarks remain clustered in nucleons. 
Why are quarks not seen to bind other than in clusters of three? How much 
excited would be the state comprising 3N freely movable quarks? In a two- 
vacuum model this is a state that should be nearly stable. But it is hard to 
imagine formation of dense quark matter without allowing for heating during 
compression in collision of heavy nuclei. At finite temperature we have some 
quark-antiquark pairs and real gluons appear. This is the quark-gluon plasma 
(QGP) state of matter last seen in the early Universe. In the limit that color 
charge interactions of quarks and gluons are small, which we believe to be 
the case in non-Abelian gauge theories at high energies, given the decrease in 
strength with energy scale of the coupling constant, the relativistic quantum 
gas properties only depend on the number of excitable degrees of freedom, that 
is the degeneracy g and the kinematic relation between energy and particle 
momentum. 

It is of course far from certain that the QGP can or/and will be created 
in relativistic nuclear collisions which last just an instant, as short as light 
needs to cross a nucleus, 10~ 22 s. Our hope and expectation is that in a 
statistical system with many degrees of freedom at as high a temperature 
as kT = 150-250 MeV, the detail of hadron dynamics and structure that 
escaped our attention will become irrelevant and we will be able to observe 
properties of deconfined nuclear matter. It is possible that the complex aspects 
of interactions within confined bound states become irrelevant in the high 
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density/energy limit. The question then is under what conditions we actually 
encounter this asymptotic limit, and what, if anything, we can learn when 
studying the approach to this limiting case. 

We shall address here several issues that we left open in our report made 
at the last meetingLl and we shall attempt to make an elementary presentation 
accessible in a large part to students; however, we recommencLthat the second 
half of this paper be read in conjunction with our earlier reportQ, as it continues 
and updates the developments reported. Before proceeding with more rigorous 
discussion of the vacuum and hadronic structure in Section ||, we shall give 
next a brief historical reminiscence about relativistic nuclear physics. Recent 
aspects of quark-gluon plasma studies are introduced in Section [| We describe 
the considerable progress that has occured in past tow years regarding strange 
particle flavor production in Section |]. The still preliminary discussion of 
particle production in hadronizing QGP is offered in Section ^, where we also 
compare with recent Pb-Pb 158 GeV A CERN heavy ion collision experimental 
results. We conclude this report with a discussion of the ongoing research work. 

Historical remarks 

The discovery of natural radioactivity more than 100 years ago occurred when 
both quantum mechanics and relativity, the two pillars of our current quest 
for quark physics, were not yet formulated in their final form. Relativity and 
nuclear physics evolved initially without mutual interaction. Quantum physics 
was applied rapidly to understand puzzles of nuclear physics, a good example 
is the a-decay using quantum tunneling. Relativity enters with Fermi theory 
of /3-decay: '...in order to obtain relativistically invariant form ... necessary 
at the velocities of emitted electrons close to velocity of light... we must use 
Dirac four functions...' writes George Gamov in Structure of Atomic Nuclei 
and Nuclear Transformations published 60 years ago. This is perhaps the first 
textbook mention of relativity within nuclear physics. 

Hideki Yukawa proposed a theoretical yet undiscovered particle, the meson, 
as the origin of short range nuclear interactions. This bold step generated a lot 
of interest in fundamental understanding of nuclear forces and nuclei, which 
continues to this day. After the initial confusion caused by the_unexpected 
heavy electron, the muon, the discovery of the 7r-meson in 1947B by Lattes, 
Occhialini and Powell just 50 years ago today made the odd couple, relativity 
and nuclear physics, inseparable. With this discovery a novel interpretation 
of the nucleon-nucleon interaction came immediately within range. Already in 
1948, Rosenfeld's book Nuclear Forces describes the Breit-type reduction of the 
fundamental meson interaction to the effective nonrelativistic form, involving 
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the spin-orbit coupling. 

Since there are several mesons, generalization of Yukawa ideas involve 
nucleon-nucleon interactions with different Lorentz transformation symmetry, 
from which the nuclear shell model potential arises, given that the range of 
the interaction is the mean distance between nucleons in nuclei. In magni- 
tude, the effective single particle nucleon potential V e g in the nucleus is not 
much bigger than a few percent of nucleon mass ton , as is born out by Bethe- 
Weizsacker mass formula which gives the bulk nuclear binding energy at about 
15 MeV/nucleon, less than 2% of ton- However, the spin-orbit coupling influ- 
ence on nucleon energies is in comparison surprisingly big, as is born out by 
study of nuclear spectra. Relativity provides here a rather simple explanation 
of a fundamental puzzle of nuclear physics, and also explains how it can happen 
that strongly interacting nucleons are bound by a relatively small potential. 
At this point relativity becomes inseparable part of nuclear physics. 

The solution of this puzzle requires consistent Breit reduction of relativistic 
wave mechanics to non-relativistic limit, schematically: 

V ce ^ { -U + V } + ^± { U + V}L.S + .... (1) 

We see that the spin-orbit force is the sum of the gradients of the (pseudo)scalar 
U and vector meson V exchange potentials, while the radial part of the effective 
potential for the nucleons is the difference of both these contributions V — U. 
This than explains the relatively small nuclear shell potential, now understood 
as the difference V — U, as U, V can be large as long as their magnitude is 
similar. Indeed, we can have complete absence of spin-orbit coupling when 
the V, U potentials add and thus the L ■ 5-coupling term vanishes, while the 
opposite limit in which the (nearly) exact cancelation of the V 1 U potentials, 
and the associated maximum of the L ■ S occurs, is the interesting property of 
nuclear interactions. Today, we believe that these nuclear interactions are the 
'unscreened', van der Waals-type forces between quarks confined to individual 
nucleons. Efforts continue to obtain QCD-quark physics based derivation of 
the residual nuclear interactions. 

The development of nuclear matter theory, neutron matter in neutron stars 
ensued in parallel with studies of high energy cosmic ray reactions. These two 
initially separate developments were precursors to our present day interests in 
hot nuclear matter, and ultimately hot hadronic matter, that is nuclear mat- 
ter at so high temperature that it contains a significant meson abundance. On 
the other side, there was early recognition among particle physicists that ele- 
mentary collisions involving strongly interacting particles lead to multi-particle 
production which could be described as if originating from a fireball of dense, 
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hot matter. Such theoretical work on multi-particle production by E. FermiB 
in USA, and L. LandauH in Russia, which paved the way to the development 
in early sixties of statistical bootstrap model description of hadron production. 
Rolf Hagedorn, working out of CERN, has been since 1964 implementing^ the- 
oretical picture of particle production from boiling hot hadronic matter E3. He 
was able to describe many experimental features of hadronic spectra within the 
newly developed statistical bootstrap model of dense hadronic matter compris- 
ing a resonance gas of point-like particles. Not only the properties of hadronic 
matter but also the hadronic mass spectrum emerged correctly from this ap- 
proach. Ultimately, a connection has been established between particle and 
nuclear hot hadronic matter descriptions, by way of introducing into Hage- 
dorn's bootstrap aside of 'elementary' mesons, also 'elementary' baryons of 
finite size til. 

In parallel to the development of phenomenological hadronic matter the- 
ory the fundamental understanding of hadronic structure was emerging. Just 
more than 25 years ago the present day stage was set when quantum-chromo- 
dynamics (QpD) was recognized as the fundamental gauge theory of strong 
interactions □ from a number of other attempts to explain hadronic structure 
in terms of quarks. The development of hadronic structure within the MIT- 
bag model followed without delay Lj, with constituents of nucleons, up and 
down quarks having in this approach just very tiny masses. To explain why 
the nucleon is so much more heavy than an electron an old model of the neu- 
tron has been resurrected: before the neutron discovery, the atomic nucleus 
was constructed from protons and tightly bound additional electrons. Without 
quarks this model could not work well since electrons were shown not to have 
strong interactions. Replacing the electrons with quarks which have indeed 
both strong and electric charge allows to describe the large masses of protons 
and neutrons in terms of the light elementary object, the 'electron' of QCD, 
the quark. 

With light quarks as building blocks of hadrons, it took no time at all 
for Garruthcrs to propose the existence of "Quarkium: a bizzare Fermi liq- 
uid" i3. Hagedorn extended his statistical bootstrap by considering the interior 
of hadrons, which were now understood as having a finite volume filled with 
light quarks. Under Hagedorn's guidance, I begun in 1977 at CERN to explore 
the melting/dissolution/fusion of dense and hot hadronic bootstrap matter 
into quark matter. One of the key issues is, if at the conditions of density and 
temperature, when confined, nucleon/meson type matter, reach the crossover 
point to the color conductive conditions of the vacuum we actually will en- 
counter a phase transition. We soon arrived at a detailed description of the 
first order phase transition that we believed occurs between the confined and 
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deconfincd phaseta. Though this subject has since many times been revisited, 
this extension of the statistical bootstrap model to finite size hadrons with 
quark internal structure is still today the foundation of all detailed models of 
the dense hadronic matter formed in relativistic nuclear collision. 



First conclusions: relativistic nuclear collisions 

Looking back to these pioneering times of quark matter days the text of the 
abstract of the inaugural lecture of June 1980, which I presented at Frank- 
furt University comes to mind. It is quite short: Quark Matter - Nuclear 
Matter. 'The fusion of the constituents - quarks — of protons and neutrons 
into quark matter, a new phase of nuclear matter, is expected to occur under 
experimentally accessible conditions of pressure and temperature.' 
Our objective is to 

- break up matter, freeing quarks, in laboratory experiments; 

- study the Universe about 40/xs after the Big-Bang by measuring properties 
of this hot, deconfmed state; 

- create a firm foundation of understanding of strong interactions in terms of 
a fundamental theoretical paradigm, confirmed experimentally. 

Nearly certainly, there is some not yet understood principle at work that is 
needed to eliminate unobserved hadronic states. This is at the same time a 
challenge and an undesirable complication within this research program. 

2 Quark Bag Model 
Vacuum Structure 

Among the most far-reaching developments of the recent 25 years of research 
into the consequences of fundamental interactions is recognition that the true 
physical vacuum is a state of considerable complex and physically significant 
structure. While we know that strictly speaking the vacuum is empty, its 
quantum structure (wave function) can be highly non-trivial, deviating con- 
siderably from that of a non-interacting Fock space. This is not the place to 
present the multitude of phenomena that go along with this effect, but we can 
illustrate some simple properties in order to justify the principal of confinement 
and the related appearance of light relativistic quarks in the description of the 
deconfmed QGP. Indeed, one finds a quite remarkable comment to this point in 
Weinberg treatiseH (Volume II, pl90, bottom) on Quantum Theory of Fields: 
'. . . this work was done without a specific theory of the strong interactions. One 
of the reasons for the rapid acceptance of quantum chromodynamics in 1973 as 
the correct theory of strong interactions was that it explained the SU(2)xSU(2) 
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symmetry [inherent in Adler-Weisberger sum rule of 1965] as a simple conse- 
quence of the smallness of the u and d quark masses.' So if quarks have a small 
mass, what is the origin of the 1 GeV scale of the mass of nucleons? 

The infrared QCD instability leads to the appearance of a finite glue 'con- 
densate' field i.e vacuum expectation value (VEV) of the gluon field-correlator 
which is evaluated from experimental data to be: 
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-G 2 ) ~ (2.3±0.3)l(T 2 GeV 4 = [390 ± 12MeV] 4 , (2) 



here a s = g 2 /47r is the strong interaction (running) coupling constant, see 
Section ||, and 

\g 2 = \ J2 G^Gr - £[B Q 2 - E 2 } , (3) 

a a 

with a — 1, . • • , N 2 — 1) gauge field components in SU c (iV). The value (and 
space-time shape) xif the glue condensate has also been extracted from lattice 
gauge calculations!!! It agrees HieLLwith the numerical value shown in Eq. (^J) , 
obtained from QCD sum-rules tBt3. Note that Lorentz and gauge symmetry 
of the vacuum state dictates that the VEV of a product of two field operators 
of the type shown here satisfy: 

(G a ^(x)G b pa (x)) = {g w9va - g ^g up )S ab ±(G 2 (x)) , (4) 

Taking the required contractions and using Eq. (||) one finds: 

(J2^) = -(J2E a 2 ) (5) 

a a 

Eq. (||) implies that (B 2 {x)) is positive, and (E 2 (x)) is negative. However, 
these signs arc defined with reference to the perturbative state \P) . Specifically, 
we are considering products of field operators which are normal-ordered with 
respect to \P) - where the condensate fields thus vanish by definition. The 
interpretation of Eq. is that the B-field fluctuates in the true QCD vacuum 
\V) with a bigger amplitude than in the perturbative 'vacuum' |P), while the 
E-field fluctuates with a smaller amplitude than in \P): \V), the true vacuum, 
may be completely 'magnetic', without any electric fluctuations, while \P) 
may be completely' electric', without any magnetic fluctuations. There are 
obviously many different ways to model and understand the glue condensate, 
and we shall not discuss this here in any greater detail. Within such models a 
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clear connection between the condensation energy B and the glue condensate 
(^f G 2 ) usually arises. 

Another important property of the true vacuum has been shown well be- 
fore the fundamental degrees of freedom and QCD have been understood. 
Expressed in modern language, the GOR-current algebra relation reads (see 
e.gl9): 

mlfl ~ --(m u + m d ){uu + dd) , (6) 

where u , d are the spinor field operators representing the two light quark flavor 
fields. The pion decay constant is / ff = 93.3 MeV, and m n refers to the charged 
pion mass m 7r + = 139.6 MeV. The non- vanishing light quark masses m u + 
are recognized to be the source of a finite mass of the pion; for vanishing quark 
masses the hadronic chiral SU(2)ixSU(2)fl symmetry is realized in QCD is 
exactly. Current .b^L estimate of the running QCD-current quark masses at 
lGeV scale are 000: 

(m„+m rf )|iGeV - 14.7 ±0.8 MeV, m s | 1GcV =s 195 ± 12 MeV. (7) 

This light quark mass estimate arises from the COR relation Eq. (|^) and the 
quark condensate deduced from the properties of the vacuum: 

i (UU + dd) |iGcV = ^ (qq) 1 1 GeV = - [(225 ± 9) MeV] 3 . (8) 

Model calculations employing mean field configurations of gauge 

fields in the QCD vacuum invariably suggest that it is the presence of the glue 
field condensate which is the driving force behind the appearance of the quark 
condensate. For example a specific modeled, (which employed a self-dual co- 
variantly constant field) for the non-perturbative gauge field configurations in 
the structured QCD vacuum, finds that the quark condensation is a minor and 
stabilizing contributor (6%) to the vacuum energy due in its bulk part to the 
glue degrees of freedom. We further note that at high temperature the vac- 
uum structure of QCD, as expressed by Eq. (|^) in. terms of the glue condensate, 
melts and one reaches the perturbative vacuumu. This confinement to decon- 
finement transformation and the chiral symmetry restoration as expressed by 
the melting of the quark condensate are seen at the same temperature. 

It is thus generally assumed that the glue condensate is driving the quark 
condensate, and thus the large chiral symmetry breaking, despite the small- 
ness of the quark masses. This is tantamount to the picture of quarks confined 
by glue fluctuations and the resulting large relativistic quark confinement en- 
ergy. We now illustrate how presence of (fluctuating) gauge fields can induce 
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appearance of Fermion condensates. Scrpiinger in his seminal paper on gauge 
invariance and vacuum fluctuations [seecj, Eq. (5.2)] shows: 



(j(xW(x)) - -mx)^{x)\) = ±J± (9) 

The left-hand-side of Eq. (^) defines here more precisely the meaning of the 
quark condensate in terms of the Fermi field operators at equal space-time 
point. The right-hand-side refers to the effective action density r[A M ] of 
Fermions in presence of gauge potentials A„. 

The one loop the effective action r' 1 ' for constant gauge fields in the 
Abelian (QED) theory was already presented by Euler, Heisenberg, KockelEH. 
We re-express this result using the invariant Maxwell-field-like-quantities E, B, 
which are related to the Maxwellian electric E and magnetic B fields by: 

B 2 = ^ (E 2 - B 2 ) 2 + 4(E ■ B) 2 - ^-(E 2 - B 2 ) (10) 
-> \eB\ 2 , for|£|^0; 

E 2 = ^(E 2 - B 2 ) 2 + 4{E ■ B) 2 + ^-{E 2 ~ B 2 ) (11) 
-> \eE\ 2 , for|B|->0. 
We also recall for later convenience: 

E 2 - B 2 = --F UU F^ U = --F 2 , E ■ B = — F UV F^ V = — FF , (12) 

2 2 4 M 4: y ' 

where: 

F"" = d^A" - d v A» , F^ = ^ vap F a0 (13) 

In this notation the effective QED action to first order in the fine-structure 
constant a e — e 2 /4ir, and evaluated in the limit that the Maxwell field is 
constant on the scale of electron's Compton wave length (which is the situation 
for all externally applied macroscopic fields) is given by: 

r(D__J_ r d _i e - 

We note the two-fold subtraction required here: the first eliminates the field 
independent, zero-point action of the perturbative vacuum, and corresponds 
to normal-ordering of the field operators. The second subtraction is absorbed 



sE 



sB 



tan sE tanh sB 



3 V 



B 2 )s 2 



(14) 
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in charge renormalization which assures that for weak-fields the perturbative 
asymptotic series begins 0(E 4 ,B 4 ,E 2 B 2 ). SchwingerES pointed out that the 
singularities along the real s-axis of the proper-time integral in Eq. ( |l4| ) are 
related to instability of the vacuum due tax pair production, a process in prin- 
ciple possible when potentials are presented that can rise more than 2m, which 
is of course the case in presence of a constant, infinite range, electrical fields. 
Employing Eq. (0) we obtain the condensate: 



l, r7 , , , , N1 , m m 2 f°° ds 
AtT 2 



m^(*),V« (1) = £2 / 



sE sB j 



tan sE tanh sB 



(15) 



We omitted here the renormalization of charge subtraction term [second suht. 
traction in Eq. (|l5|)1 needed only for charge renormalization. Schwinger cj 
reintroduced this specific term to obtain the leading effective electromagnetic 
interaction coupling of a neutral scalar meson to two photons. 

There are few if any observable macroscopic effects of the fermion vacuum 
fluctuation in QED as these need to be induced by an extremely strong external 
force. In QCD situation is fundamentally different, due to spontaneous field 
fluctuations in the glue sector. To obtain a measure of the effects we need to 
substitute in the effective action Eq. ( pT| ) the pure U(l) Maxwell gauge field 
by the invariant combination of U(l) and the condensed SU(3) gauge fields: 

^ { ^ G 2 } + ^ q 2p2 (16) 
TT TT TT 

Here q — 2/3, —1/3 are the quark fractional Maxwell charges. In general one 
can 'forget' about the QED charges and proceed with QCD fields G alone. It 
is not appropriate to expand the effective action perturbatively since the QCD 
fields, even in the vacuum condensate, are in general more significant than the 



masses of light and strange quarks. When studying Eq. (15) one is rapidly lead 
to the conclusion that the mass of light quarks is not the important scale in the 
problem at hand, since for u, d, s quarks the glue condensate scale is greater 
than the quark mass. 

We note that the substitution we have made in Eq. ([n|) implies when 
effected in the highly nonlinear function Eq. (|l5|) that we have replaced 

(G 2n ) -» ((G 2 ))™ (17) 

which is correct under the tacit assumption that the fluctuations of the glue 
con den s ate C an be understood as if generated by a background stochastic 
fieldEZroEl This is than the picture of the true QCD vacuum that we consider 
to be very promising model of the complex reality and which is implied in the 
following. The stochastic field fluctuations confine the color charges of quarks 
and gluons, and determine the confinement size. 
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Hadronic (Hyperfine) Structure 

Without going much into well studied detail, the simple yet successful picture 
of hadronic structure is the 'bag'-model. We view a hadron as comprising a 
spherical confining volume of radius R containing iV-quarks in a bound state 
of mass Mn- Neglecting here the quark mass compared to the zero-point 1/R 
energy the bound state mass is given by: 

M N = — + -kR 3 B + Hi . (18) 
K o 

The residual QCD interaction term Hi comprises the quark-quark interac- 
tions. The energy eigenvalue n/R of the relativistic quark wave function with 
a boundary condition enforcing a vanishing current through the surface gives 
K = K c = 2.04. The second term in Eq. ( |T8|) is the bag volume energy B, the 
latent (condensation) heat of the vacuum. The balance of forces requires that 
the stationary state is at the minimum of Mff(R) with respect to R, whence 
we determine the mass of the V-quark bound state: 

M N = ( K iV) 3 / 4 1.755£ 1/4 for Hj = 0. (19) 

If we fit the nucleon mass M 3 = 940 MeV by Eq. @, than we find B —> 
B^ 4 = 140 MeV, but working out the A- mass, M 3 = 1330 MeV one finds 

1 /4 

S A ' = 190 MeV. These two values differ by factor 3.4 and thus we have to 
refine the approach. What is missing is the strong quark-quark interaction that 
is capable to generate the hadronic multiplet mass splitting. The color-electric 
quark-quark interaction is comparatively weak as the overall color charge of 
all hadronic states is zero, and only particles comprising both light and heavy 
quarks have a sizable static electric color charge density distribution. The 
color-magnetic hyperfine 'spin-spin' interaction term isll3: 

Hj = J2 Y^"^» (20) 

i>j£h 

Here S — er a /2, a = 1,2,3 is the spin vector of the two interacting quarks 
i, j, written in terms of the Pauli spin matrices a a , and A a /2, a = 1 . . . 8 is the 
generator of the color non-Abelian SU(3) charge. The dot product indicates 
summation over a = 1,2,3 for SU(2) and a = 1 ... 8 for SU(3). 

It turns out that all strongly interacting particles known (mesons and 
baryons such as n, p, N, A, K, A)j-can be consistently fitted using this color- 
magnetic interaction Eq. ( pp| ) , see Let us recall here a ab-initio fit to all 
strange and non-strange mesons and baryons, allowing for massive strange 
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quark variable scaled which produces B^/ 4 = 170 MeV. In this fit also the 
quark eigenenergy k was fitted, and the fitted value k — 1.97 is nearly the 
theoretical value nth = 2.04, the small difference between theory and fit ex- 
pected to arise from the distortion of the wave-function due to the hyperfine 
color-magnetic interaction. The spectra of hadrons tell us that quarks live in a 
relativistic, confined orbital, which fact supports the QCD-quark-confinement- 
vacuum structure picture of the nucleon. 

Another success of both the relativistic and the non-relativistic quark bag 
model of hadrons is that we find that the most stable bound quark states 
have quantum numbers seen in stable hadrons, and that the mass hierarchy 
of flavor mass multiplets is correctly explained. This provides firm evidence 
that the interaction between quarks is obeying the rules of non-Abelian color- 
SU(3) algebra. To see this let us consider in more detail the situation for three 
quark states, i.e. baryons. We compute the strength of the chromo-magnetic 



interaction fi 1 } in Eq. (20), which we define as follows: 

H I = 8E l i h I ; mN- E (My-y^'K)- (21) 

Signs are chosen such that 5E cx 1/R is a positive quantity, fij can be positive 
or negative, depending on the quantum numbers of the bound state. 

To obtain the value of the matrix element uli in a given hadronic state we 
follow here the permutation operator method E3. The exchange operator Pij 
of two quarks is composed of three factors for spin, color and flavor quantum 
numbers, assuming that in the ground state the spatial wave function is iden- 
tical for all quarks. For the totally antisymmetric quark state \h q ), for each 
quark pair i,j we have thus 

P?jPtj\h) = -P&h); p Z = l + 2 Y'Y ; Pti = l + *S i >S* (22) 

The above explicit form of the permutation operators for the SU(2)-spin and 
SU(3)-color groups follows from the commutation properties of the generators 
S = <r/2,A/2. Inserting Eq. (^2|) into Eq. (pl|) we obtain for the hyperfine 
interaction operator: 



f 1 , 1 1 X 1 X j S i ■ S j I 



(23) 



and naturally /ij = (fij) . 

To complete the evaluation of fii we need to commit ourselves to the 
symmetry group governing the flavor exchange operator. Since strangeness is 
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well distinguishable from u, d flavors one can take the point of view that we 
should use only the light quark SU(2) flavor group with: 

Pi- = J + 2/ 1 • P , (24) 



2 



and hcnccli 



P 12 ^ \ 2 6 4 2 2 J ' v; 

where we have used Yli>j = n ( n ~ I)/ 2 - We than have to account for the 
interaction of strange quarks with the non-strange quarks separately. It turns 
out that the fittedEj chromo-magnetic interaction energy SE qs between strange 
and light quarks is half as large compared to u, d interactions, and that strange- 
strange chromo-magnetic interaction 5E SS is yet five times weaker than the u, d 
case. Such decrease of magnetic interaction is consistent with the expectations 
that the more massive quark has a smaller magnetic moment. While these 
results support the need to separately treat the light and strange flavors, we 
note in passing that if we had taken u, d, s quarks at par, than the larger SU(3)- 
flavor group would relate to flavor symmetry. Accordingly in that case one 
should evaluateE-3 in Eq. ( 23 ) the flavor exchange operator using the quadratic 



Casimir of the SU(3)-flavor group to represent the flavor permutation operator 
in Eq. (p2|). We will not pursue this approach further here, which has merit 
when considering other aspects of hadronic spectra. 

Using the SU(2) flavor exchange operator we find that the chromo-magnetic 
interaction term of light quarks u, d assumes the form: 

A 1 = "' ( "V 6) + 3^ +v + h Sl{Sl + + q) ■ (26) 

where the index Zight reminds us that here we have only accounted for the 
light quark u, d contribution and quantum numbers in all hadrons. In deriving 
Eq. fl2q ) we have used some well known relations: 



i>j i j i 

2^.^" = (]r^).(E^)-E^: 



i>j 



2 E A4 - aj = (E A4 hE aj )-E ai - ai 
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and 

C£ n 2 \h) = Hi + i)\h) , (£ si )» = s ( s + ■ 

i i 

given the total spin and isospin operators introduced implicitly above. Ev- 
erything works in the same way for the color group SU(3), except that the 
individual quadratic Casimir operator has a different eigenvalue to keep in 
mind. In the fundamental representation 




which is just the inverse of 3/4, the eigenvalue of the quadratic Casimir opera- 
tor for the SU(2) in the fundamental representation 2. For the composite state 
\h) of particles q and antiparticles q in fundamental representation we have: 

E(Et) \h)=C 2 (p,q)\h), 

a=l \ i / 

where the values of (p, q) for a few representations are: 1=(0,0), 3=(0,1), 
3=(1,0), 6=(2,0), 6=(0,2), 8=(1,1), 1Q=(3,0), etc. More generally, the dimen- 
sion of the representation and the Young-tableaux numbers p, q are related 
by: 

d(p,q) = i(p+l)(g + l)(p + g + 2). 
The quadratic Casimir operator in terms of p, q assumes the values 

C2(P, 9)=P + 9+ y + y + y- 

Explicitly: 

C 2 (1) = 0, C 2 (3) = C 2 (3) = ^, C 2 (6) =C 2 (6) = 3^, 0,(8) = 3. 

It is with great satisfactions that we note that Eq. ( |2rj ) has a minimum 
for ni = 3, corresponding to the baryon system, and that overall the chromo- 
magnetic interaction energy factor is negative only for the nucleon state among 
the three quark states. The splitting of TV and A (difference between I = 
3/2, S = 3/2 and / = 1/2,5 = 1/2 states) M A - M N = 6E thus leads tog 

a When fitting the mass of the A one has to note that the coupling to a strong decay channel 
also shifts the quark based hadron mass downward yielding the physical mass, and hence 
the quark-based portion of the A mass is to be taken at + Y & 
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SE ~ -400MeV, which is consistent with the splitting M s - M A = (3/16) SE 
(difference between I — and I — 1 states). 

These are very persuasive examples of the successes of the bag model of 
quark structure of hadrons, and prove that the hadron masses are properly 
described by the color magnetic interaction. 

Limitations of the quark-bag picture 

But if Eq. (|2^) works very well for the fundamental multiplets of baryons and 
mesons, why is it that it also predicts unwanted states such as baryonium qq- 
qq, the qqs-qqs di-baryons and many more, not seen? Some of these exotica 
are quite stable: which multi-quark states are most hmmd can be deduced 
from the general bonding (Hund) rule of quark physics E3: 

• the quarks and antiquarks are separately in the largest possible repre- 
sentation of color and spin, and 

• the total state is in the smallest possible representation of spin an color. 

Thus the state with the structure [(uds)s r A s 8> (wds)8 ( .,4 s ]o c ,o s where the sub- 
indices refer to color and spin multiplicities of the representations is the most 
bound exotic 6 quark state. Such a hadron has not been found despite literally 
twenty different experimental searches. So probably not all is well with our 
current understanding of quark structure in hadrons and strong interactions. 
It could be that the small size of normal hadrons may not yet require a fully 
dissolved vacuum structure, while the larger exotica do require full structure 
dissolution, In such a case there is additional vacuum energy cost required 
to form these larger exotic states, a fact which would induce considerable 
instability by hadronic dissociation. The hope is that the study of really large 
chunks of quark matter we call quark-gluon plasma will provide a shortcut to 
a resolution of this problem. 

_Lct us hence make a step back and look again at the hadron structure 
fit E3 - if this fit is so successful, what does it tell us about the volume energy, 
which is the 'condensation' latent energy of the vacuum? The fitted value 
Bl /A = 170 MeV, which we also can express as 2?b = 0.1Gev/fm 3 is relatively 
small. The magnitude derives from the fact that the quark-related energy 
component in the hadron is 3/4 of the total, as can be seen solely on dimen- 
sional grounds inspecting Eqs. (|l^;|l^), since only terms changing with R as R 3 
and 1/R appear. The volume energy component in the nucleon is thus about 
to/4 = 235 = BhVh MeV. Only if we could make hadron volume Vt very small, 
the volume energy could come out big. However the elastic electromagnetic 
form factors and other electromagnetic properties of hadrons indicate that the 



1G 



Maxwell-charge distribution inside hadron is not less than 0.7-0.8 fm in size. 
The above considered value of the volume energy leads to a quark confinement 
volume with radius R — 0.83 fm, for protons, compared with the charge radius 
< r 2 >= (0.86 fm) 2 . Clearly, there is not much room for change and surely we 
will not be able to find B = 0(1) GeV/fm 3 , a value motivated by numerical 
simulations of lattice QCDH. 

3 Quark-Gluon Plasma 
Relativistic heavy ion experiments 

In order to form relatively large quark-gluon filled space-time regions the best 
tool we have today are large nuclei. These nuclei are made to collide at very 
high energies, many times higher than the rest mass of particles involved. The 
relativistic energy is required to produce regions of space filled with movable 
color charges of quarks and gluons, the quark-gluon plasma (QGP). Primarily 
because the identification of this new, locally deconfined form of matter is 
difficult, the object we search for seems to be the never-to-be-found nuclear 
Holy Grail. However, this deconfined QGP state must in principle exist if 
QCD is the true theory of strong interactions: such a free quark-gluon phase 
is clearly seen in numerical studies of QCD within the lattice gauge theory 
(LGT) simulations, though the details of the quark structure of the deconfined 
state remain today somewhat obscured by the difficulties that one encounters 
in numerical treatment of FermionstX 

The study of highly excited and dense hadronic matter by means of ultra- 
relativistic nuclear collisions is a relatively novel, interdisciplinary area of re- 
search in rapid experimental and theoretical evolution. It is closely related 
to the fields of nuclear and particle physics and, accordingly, our material en- 
compasses aspects of both these wide research areas. Among 'consumers' of 
the material presented here we also encounter researchers in the fields of as- 
trophysics and cosmology. The idea to test this is to squeeze and compress 
large nuclei, such that individual nucleons are made to fuse into a new common 
structure comprising freely moving quarks. 

The initial experimental programs were launched nearly 30 years ago at 
the Lawrence Berkeley Laboratory (LBL) at Berkeley, USA, and at the Dubna 
Laboratory (Russia) for lighter ions. Another 1-2 GeV facility, SIS, of com- 
parable physics perspective has come on line in recent years at GSI in Darm- 
stadt, Germany. A greater beam intensity and new detector technologies at 
SIS allow to explore systematically phenomena that could not be studied at 
the BEVALAC. Studies of properties of excited nuclear matter performed at 
these facilities display particle spectra 'temperatures' of the order of 100 MeV, 
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accompanied by relatively abundant pion production. 

These first advances where made with relatively modest energies, and the 
results had mainly bearing on the properties of the nuclear matter close to the 
conditions encountered in supernova explosions. However, they demonstrated 
the possibility to study the properties of compressed and excited nuclear mat- 
ter in laboratory. Together with the theoretical work, these results implied 
that QGP could be in reach of existing accelerator facilities. Thus a search for 
the point of transition from the hadronic gas phase of locally confined nuclcons, 
mesons, etc, to the deconfincd QGP phase that gave birth to the development 
of the present day research program, in USA at BNL (Brookhaven National 
Laboratory, Upton - Long Island, New York), and in Europe at CERN (Euro- 
pean Nuclear Physics Research Laboratory, Geneva). The first Oxygen beam 
at 60 A GeV was extracted from the SPS (Super Proton Synchrotron) ac- 
celerator at CERN and met target in the late fall of 1986, about the same 
time as BNL started its experimental program at AGS (Alternate Gradient 
Synchrotron) accelerator with 15 A GeV Silicium ions. Very soon after, the 
energy of the SPS beam could be increased to 200 A GeV and a Sulphur ion 
source was added. In order to allow experiments to study the relatively large 
volumes and longer lifetimes expected in dense matter formed in collision of 
heaviest nuclei, an upgrade of SPS injection system was carried out, which 
allowed to accelerate the Lead ions to 158 A GeV in the Fall of 1994. At BNL, 
a similar development allowed to accelerate Gold beams to 11 A GeV at the 
AGS. 

Today we are repositioning our interests around two new experimental 
facilities: RHIC and LHC. The RHIC (Relativistic Heavy Ion Collider) at BNL, 
is slated for physics progress in 1999, with colliding nuclear beams at up to 100 
A GeV, thus allowing exploration of a entirely new domain of energy, 10 times 
greater than currently available at SPS. The LHC (Large Hadron Collider) 
which is being developed at CERN, should be commissioned at midpoint of 
the first decade of XXI century, and it will devote an important part of its 
beam time to the acceleration of nuclear beams. The final energy available at 
the LHC is expected to be a factor 30 or more higher compared to RHIC. 

Many experiments arc in progress both at CERN and BNL and our present 
interest addresses primarily results obtained at the higher CERN energies, 
involving strange particle production. Several experiments comprise similar 
'crews' and we will group their names together. This leaves us with two major 
lines of approach: NA35/NA49 and WA85/WA94/WA97/NA57. While the 
NA35 etc groups are interested in obtaining full phase space coverage and use 
detectors such as streamer chambers and time projection chambers, the WA85 
etc groups concentrate on a small window of phase space 'opportunity' with 
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a spectrometer, which allows to identify and observe rarely produced multi- 
strange particles. Another series of experiments which is of some importance 
for our study of heavy flavor quarks is NA38/NA50/NA51 series, which ob- 
serves by means of their dimuon decays the production of vector mesons. 

The first conclusion one can draw from the recent results including the 
collisions of Pb-PB nuclei involving projectiles colliding with a laboratory fixed 
target at 158 A GeV is that there is a localized space-time region in which 
high concentration of matter and energy is reached: we refer to this object 
as 'fireball'. We believe that local thermal equilibrium is established within 
the volume occupied by the fireball, and thus particles emerge with spectra 
characteristic of the surface temperature and velocity. 

Probes of dense hadronic matter 

Since in the collision of large nuclei, the highly dense state is formed for a 
rather short time of magnitude 2R/c, where R is the nuclear radius, one of 
the major challenges has been to identify suitable physical observable of de- 
confincment. This difficult problem of detecting reliably the formation of an 
unknown phase of matter, existing only as short as 0.5 • 10~ 22 sec, has not 
been completely resolved today. The electromagnetic probes involving directly 
produced photons and dileptons are witnesses to the earliest moments of the 
reaction, but their production rates are in general very small, the experimental 
7 yield is dominated by the secondary processes 7r° — * 7 + 7. The dilepton 
yield is also mostly resulting from meson decays, but there are small kinematic 
regions where this background yield is very small. Moreover, the dilepton 
spectrum gives interesting insights about the vector meson yields and their 
variation with experimental conditions. We will briefly address here this very 
interesting observable. 

The principal observable of the relativistic nuclear collision process are 
spectra of hadronic (strongly interacting) particles. Their multiplicity is grow- 
ing rapidly with both energy and size of the projectile-target nuclei. Aside of 
the single particle spectra one can also relatively easily measure two particle 
correlations. The HBT correlation measurement is widely used to determine 
the geometric properties of the central fireball, and the results agree with a 
reaction picture between the nuclei relying on geometric considerations. How- 
ever, this approach is still plagued by fundamental and unresolved issues. Work 
is progressing rapidhiin this field and which is addressed in full in this volume 
by Baym and Heinz E3. 

Strangeness (and at RH1C/LHC also charm) and entropy c3 are good 
hadronic observable since both will be preserved by 'reasonable' evolution sce- 
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narios of the dense matter fireball: the melted QGP state is in general entropy 
richer than the frozen HG phase. Once entropy has been generated, it cannot 
be lost, an entropy excess accompanies QGP formation. When abundant par- 
ticle production is possible, this entropy excess is seen as an enhancement in 
the total hadronic particle multiplicity with each (relativistic) meson carrying 
about 4 units of entropy out of the interaction region. Similarly, strangeness is 
in general more abundant in QGP than HG phaseELand it is not re-annihilated 
in rapid decomposition of the dense matter state Ej. It has become a key di- 
agnostic tool of dense hadronic matter because: 

1) particles containing strangeness are produced more abundantly in rela- 
tivistic nuclear collisions than it could be expected based on simple scaling of 
p-p reactions; 

2) all strange quarks have to be made, while light u, d quarks are also 
brought into the reaction by the colliding nuclei; 

3) because there are many different strange particles, we have a very rich 
field of observable with which it is possible to explore diverse properties of the 
source; 

4) theoretical calculations suggest that glue-glue collisions in the QGP 
provide a sufficiently fast and thus by far, a unique mechanism leading to an 
explanation of strangeness enhancement. 

There are two generic flavor observable (strangeness and charm) which we 
study analyzing experimental data, and we introduce these here, without an 
effort to 'orthogonalizcd', i.e., make them independent of each other: 

• absolute yield of strangeness/charm 

Once produced in hot and dense hadronic matter, e.g., the QGP phase, 
strangeness/charm is not rc- annihilated in the evolution of the decon- 
fined state towards freeze-out, because in the expansion and/or cooling 
process the rate of production/annihilation rapidly diminishes and be- 
comes negligible. Therefore the flavor yield is characteristic of the initial, 
most extreme conditions. 

• phase space occupancy 7$ 

7i describes how close the flavor yield per unit of volume (i — u, d, s, and 
in some cases charm i = c occupancy will be also considered) comes to 
the chemical equilibrium expected; 7, impacts strongly the distribution 
of flavor among final state hadronic particles. 

Because of the high density of the QGP phase, the phase space occupancy 7, 
can saturate rapidly, and thus particle abundances will emerge from a chem- 
ically equilibrated u, d, s phase, which is hardly imaginable for conventional 
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reaction mechanisms. Because entropy and strangeness are enhanced in a sim- 
ilar way in QGP, the specific yield of strangeness per particle produced is not 
a good quantity to use when searching for the deconfined state. Many other 
strategies are available, of which we favor measurement of the specific entropy 
yield per participating baryon, accompanied by a study of relative strange an- 
tibaryon yields, involving particle ratios such as A/pEJ. It is remarkable that 
the pertinent results obtained for S-Pb collisions by the NA35 collaborations 
have shown the QGP related enhancement. 



Recent results on strangeness 

The possibility that strange particle anomalies seen in recent years at SPS in 
Sulphur induced reactions on heavy nuclei are arising in consequence to the 
formation of a deconfined QGP phase has stimulated the intense continuation 
of the experimental research program in the considerably more difficult, high 
particle multiplicity environment arising in Pb induced reactions, which are 
presently possible at 158A GeV. 

The experiment WA97 u3 has further reported several specific strange 
baryon and antibaryon ratios from Pb-Pb collisions at 158 A GeV, compris- 
ing 30% of inelastic interactions. All ratios are obtained in an overlapping 
kinematic window corresponding effectively to transverse momentum p± > 0.7 
GeV, within the central rapidity region y <= y cm ± 0.5. They have been cor- 
rected for weak interactions cascading decays. The experimental values are: 

Ra = j =0.14 ±0.03, i? s = ^ = 0.27 ±0.0, Ra= ^ = 0.42 ± 0.12, (27) 

i?P = v = 0.14 ±0.02, R? = = = 0.26 ± 0.05 , (28) 

A A 

K P = § = 0.19 ±0.04, R s p = = = 0.30 ± 0.09. (29) 

Here, the lower index s, resp. s, reminds us that the ratio measures the density 
of strange, resp. anti-strange, quarks relatively to light quarks. The upper 
index p indicates that the ratio is taken within a common interval of transverse 
momenta (and not common transverse mass). We compare Rg results with 
earlier measurements in Fig.[j]. The strange antibaryon enhancement effect is 
re-confirmed in the Pb-Pb data, and we see that there is no major change 
of this result, which determines the phase space occupancy of strangeness, as 
we move from S-S or S-W/Pb results to Pb-Pb results. There is agreement 
between WA97 and NA49 on the value of Ra, even though the data sample 
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Figure 1: Sample of World results (as 'function' of experiment name for yields involving 
ratio of strange to non-strange quarks in baryons. Dark squares: recent Pb-Pb results. 



of NA49 is taken for more central trigger, constrained to as few as 4% of 
most central collisions. The—cuts in p± and y are nearly identical in both 
experiments. From Fig. 3 in 113, we obtain the value R\ = 0.17 ± 0.03, which 
we shall combine with the value given by WA97 and we thus take in out data 
fit: 



Ra = j =0.155 ±0.04. 



(30) 
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4 Strangeness Production 

The coupling constant of strong interactions 

The observed strange particles are presumably born in a deconfined phase 
which can make strange quark pairs effectively. It turns out that the production 
of strangeness is originating primarily in gluon-gluon collisions. We use two 
particle collision processes to evaluate thermal flavor production in QGP, as 
described in a. However, since this report we have realized c3e3 that recent 
precise measurements of the strong interaction coupling constant a s allow us 
to eliminate a lot of arbitrariness from the earlier calculations regarding the 
strength of the collision cross sections. The experimental measurements of a s 
occurred at the scale of the the Z mass, thus 92 GeV, and hence one has to use 
the powers of the renormalization group of QCD to evolve these results towards 
the here relevant domain. In the process we also allow for the evolution of the 
strange quark mass, which is also significant and is driven by the evolution of 
the coupling constant. We will now briefly explain how this is accomplished. 

To determine the two QCD parameters required, we will use the renormal- 
ization group functions f3 and 7 m : 

= P(<Xb(p)) , Mg- = -mj m (a s (fi)) . (31) 
For our present study we will use the perturbative power expansion in a s : 

/3 pcrt ~a 2 s [b + b ias + ...], 7 P° rt = a s [ c + Cl a s + ...], (32) 

For the SU(3)-gauge theory with ri{ Fermions the first two terms (two 'loop' 
order) are renormalization scheme independent, the three loop level 62 is known 
in both MS (minimal subtraction) and MS (modified minimal subtraction) 
renormalization schemes c3o, while C2, so far we could ascertain, was only 
derived in the MS scheme. In any case we will restrict our study to the two 
loop order which does not require renormalization scheme improvements of the 
form of the two body cross sections. To this order we have: 

60 = i( n "H' 61 = i( 51 -y" f )' (33) 
eo=f ci =iM 101 -r if )- (34) 

The number rif of Fermions that can be excited, depends on the energy 
scale [i. We have implemented this using the exact phase space form appro- 
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priate for the terms linear in rif 




n f (p) = 2+ > \n---r 1 + — e(M-2mi), (35) 



with m s = 0.16 GeV, m c = 1.5 GeV, nib = 4.8 GeV. We checked that there 
is very minimal impact of the running of the masses in Eq. (^5|) on the final 
result, and will therefore not introduce that 'feed-back' effect into our current 
discussion. The largest influence on our solutions comes here from the bottom 
mass, since any error made at about 5 GeV is amplified most. However, we 
find that this results in a scarcely visible change even when the mass is changed 
by 10% and thus one can conclude that the exact values of the masses and the 
nature of flavor threshold is at present of minor importance in our study. We 
show the result of numerical integration for a s in the top portion of Fig. |[ 



First equation in (31) is numerically integrated beginning with axuinitial value 
of a s (Mz). We use in this report the August 1996 World averageci a s (Mz) — 
0.118 for which the estimated error is ±0.003 . This value is sufficiently precise 
to reduce this uncertainty below that has limited our earlier study e9. In 
addition, the thin solid lines present results for a s (Mz) = 0.115 till recently 
the preferred result in some analysis, especially those at lower energy scale. 
As seen in Fig.^, the variation of a s with the energy scale is substantial, and 
in particular we note the rapid change at and below fi = 1 GeV, where the 
strange quark flavor formation occurs in hot QGP phase formed in present day 
experiments at 160-200 A GeV (SPS-CERN). Clearly, use of constant value of 
a s is hardly justified, and the first order approximation often used: 

_ 26 ~ 1 (n f ) 

^ = InWAoO*)) 2 ' (36) 

can be seen to be not a good approximation till rather high /i = 50 GeV scale 
is reached: we insert the numerically computed value of a s into Eq. ( |36| ) and 
obtain Ao(p) which would yield the given a s (/i). This procedure leads to a 
strongly scale dependent Aq (fi) shown in the middle section of Fig. ^. 

Given a s (u), we can integrate the running of the quark masses, the second 
equation in (J3l|). Because the running mass equation is linear in m, it is 
possible to determine the universal quark mass scale factor 

m r = m{n)/m{no) ■ (37) 

Since a s refers to the scale of [Iq = Mz, it is a convenient reference point also for 
quark masses. As seen in the bottom portion of Fig.0, the change in the quark 
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Figure 2: a s (n) (top section); the equivalent parameter Ao (middle section) and 
m r (/i) = m(jj,)/m(Mz) (bottom section) as function of energy scale /i. Initial 
value a B (Mz) = 0.118 (thick solid lines) and a s (Mz) = 0.115 (thin solid lines). 
In lower section the dots indicate the strangeness pair production thresholds for 
m B (Mz) = 90 MeV, while crosses indicate charm pair production thresholds for 
m c (M z ) = 700 MeV. 



mass factor is highly relevant in the study of strangeness and charm production, 
since it is driven by the rapidly changing a s near to ji ~ 1 GeV. For each of 
the different functional dependencies ot s (pi) we obtain a different function m r . 
The significance of the running of the charmed quark mass cannot be stressed 
enough, especially for thermal charm production occurring in foreseeable future 
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experiments well below threshold, which makes the use the exact value of m c 
necessary. 

Given these results, we find that for a s — 0.118 and m s (Mz) = 90 ± 
18 MeV a low energy strange quark mass m s (lGeV) ~ 200 ± 40 MeV, in 
the middle of the standard range 100 < m s (lGeV) < 300 MeV. Similarly we 
consider m c (Mz) = 700 ± 50 MeV, for which value we find the low energy 
mass m c (l GeV) ~ 1550 ± 110 MeV, at the upper (conservative for particle 
production yield) end of the standard range 1 < m c (lGeV) < 1.6 GeV. The 
energy threshold for pair production E\ , i = s, c, has to be determined in a 
separate calculation since it is related to the scale of energy required for the 
production of two zero momentum particles E t h — 2m(/i = -Eth), and thus: 



This effect stabilizes strangeness production cross section in the infrared: below 
-y/s = 1 GeV the strange quark mass increases rapidly and the threshold mass 
is considerably greater than m s (l GeV). We obtain the threshold values E* h = 
605 MeV for a s (M z ) = 0.118 and B B th = 560 MeV for a s (M z ) = 0.115. Both 
values are indicated by the black dots in Fig.|^. For charm, the running mass 
effect has the opposite effect: since the mass of charmed quarks is listed in 
tables for /i = 1 GeV, but the here explored value of the mass is above 1 
GeV, the charm production threshold turns out to be smaller than expected. 
There is little uncertainty with regarding the coupling constant; for m c (Mz) = 
700 MeV the production threshold is found at ~ 2.3 GeV rather than 
2 • 1.55 = 3.1 GeV, the value we would have obtained using m c (lGeV) in 
this case. This reduction in threshold energy enhances strongly the thermal 
production of charm at 'low' temperatures T ~ 250 MeV, where it is suppressed 
exponentially by a factor cx exp(— E th /T). 

QCD strangeness production processes 

The generic angle averaged two particle cross sections for (heavy) flavor pro- 
duction processes g + g — > / + / and q + q — ► / + / , are among classic results 
of QCD and are given by 



Ef 1 = 2mi{M z )m T {Ef) . 



(38) 
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where W(s) = y/l — Amffs , and both the QCD coupling constant a s and 
flavor quark mass mt will be in this work the running QCD parameters. In 
this way a large number of even-a s diagrams contributing to flavor production 
is accounted for. 

Other processes in which at least one additional gluon is present are not 
within the present calculational scheme. While only in very high density envi- 
ronment we could imagine relevant contributions from three body initial state 
collisions, emission of one or even several soft gluons in the final state could be 
relevant, thus this subject area will be surely revisited in the future. We note 
that a process in which a massive 'gluon', that is a quasi-particle with quantum 
numbers of a gluon, decays into a strange quark pair, is partially included in 
our work. At the present time we do not see a systematic way to incorporate 
any residue of this and other effects, originating in matter surrounding the 
microscopic processes, as work leading to understanding of renormalization 
group equations in matter (that is renormalization group at finite temperature 
and/or chemical potential) is still in progress e3. 

The master equation for flavor production, allowing for the detailed bal- 
ance reactions, thus re-annihilation of flavor, is: 

d -M)- = V{t)A i [l- 1 l{t)] • (41) 

where A4 is the invariant rate per unit volume and time, that is the thermal 
average of the production cross section: 

A-i — Agg + A U u + A^ + . . . 

<PP2 



1 



7;9gfg(Pl)fg(P2)crgg{s) + U { g Jq(pi)fq{p2)(Tqq{s) 



(42) 



The dots indicate that other mechanisms may contribute to flavor production. 
The particle distributions / are in our case thermal Bose/Fermi functions, and 
<7q — 6, <? g = 16 . For strangeness production nf = 2, and for charm production 
rif = 3 , will be used. 

We introduce also the flavor production relaxation time constant: 

,_i prim) (43) 

2 (Agg + A qq + . . .) 

Here rhi is the mass at the scale of energy under study; we recall that the 
equilibrium distribution is result of Boltzmann equation description of two 
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body collisions. Thus the mass arising in the equilibrium density p°° in Eq. (43) 
is to be taken at the energy scale of the average two parton collision. We adopt 
for strangeness a fixed value m s = 200 MeV, and also rh c = 1, 500 MeV, and 
observe that in the range of temperatures here considered the precise value 
of the strange quark mass is insignificant, since the quark density is primarily 
governed by the T 3 term in this limit, with finite mass correction being 0(10%). 
The situation is less clear for charm relaxation, since the running of the mass 
should have a significant impact. 




2 i i i i i i i i i i i i i i i i i 

200 250 300 350 

T [MeV] 

Figure 3: QGP strangeness relaxation time, for a s (Mz) = 0.118, (thick line) and = 
0.115 (thin line); m s (Mz) = 90 MeV. Hatched areas: effect of variation of strange 
quark mass by 20%. Dotted: comparison results for fixed a s = 0.5 and m s = 200 
MeV. 



We can now proceed to evaluate the relaxation times for strangeness and 
charm production. The current calculation offers an upper limit on the actual 
relaxation time, which may still be smaller, considering that other mechanisms 
may contribute to flavor production, as is shown by the dots in Eq. (E2J). We 
show in Fig.|| also the impact of a 20% uncertainty in m s (Mz), indicated 
by the hatched areas. This uncertainty is today much larger compared to 
the uncertainty that arises from the recently improved precision of the strong 
coupling constant determination c3. We note that the calculations made at 
fixed values a s = 0.5 and m s = 200 MeV (dotted line in Fig.[|) are well within 
the band of values related to the uncertainty in the strange quark mass. 

Since charm is somewhat more massive compared to strangeness, there 
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Figure 4: Solid lines: thermal charm relaxation constant in QGP, calculated for 
running a s (M z ) = 0.115; 0.118, (indistinguishable), m c (M z ) = 700 MeV. Lower 
dotted line: for fixed m c = 1.1 GeV, a s — 0.35; upper doted line: for fixed m c = 1.5 
GeV, a E = 0.4 . Hatched area: effect of variation m c (M z ) = 700 ± 50 MeV. 



is still less uncertainty arising in the extrapolation of the coupling constant. 
Also the systematic uncertainty related to the soft gluons (odd-a s ) terms are 
smaller, and thus the relaxation times r c we show in Fig.|4| are considerably 
better defined compared to r s . There is also less relative uncertainty in the 
value of charm mass. We also show in Fig.|] (dotted lines) the fixed m c , a s 
results with parameters selected to border high and low T limits of the results 
presented. It is difficult to find a good comparative behavior of r c using just 
one set of m c and a s . This may be attributed to the importance of the mass 
of the charmed quarks, considering that the threshold for charm production is 
well above the average thermal collision energy, which results in emphasis of 
the effect of running charm mass. In the high T-limit the choice (upper doted 
line in Fig.||) m c — 1.5 GeV, a s = 0.4 is appropriate, while to follow the result 
at small T (lower doted line in Fig.||) we take a much smaller mass m c = 1.1 
GeV, a s = 0.35 . 
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Temporal evolution of the fireball 

In order to compute the production and evolution of strangeness (and charm) 
flavor a more specific picture of the temporal evolution of dense matter is 
needed. In a simple, qualitative description, we assume that the hot, dense 
matter is homogeneous. We consider that, in Pb-Pb collisions at SPS, the 
radial expansion is the dominant factor for the evolution of the fireball proper- 
ties such as temperature/energy density and lifetime of the QGP phase. The 
expansion dynamics follows from two assumptions: 

• the (radial) expansion is entropy conserving, thus the volume and tempera- 
ture satisfy 

V ■ T 3 = Const. (44) 

• the surface flow velocity is given by the sound velocity in a relativistic gas 

«f = l/V3. (45) 

It is well understood that expansion into 'vacuum' of relativistic matter can 
proceed at any velocity, up to velocity of light. However, here Vf is understood 
not as the velocity of the first flow edge, it is the velocity of the flow of matter 
from interior of the fireball to the near surface region, which cannot exceed the 
speed of sound in the matter, unless a shock is created. We will return to this 
intricate issue in the near futureEJ. 

These two assumption imply the following explicit forms for the radius of 
the fireball and its average temperature: 

»-*■+>-*•>• r = 1+( .-IW - (46) 

We shall see below that if QGP formation is involved, a fit of strange an- 
tibaryons data either leads to direct emission before expansion, or to emission 
from a surface expanding with just this velocity V{. 

The initial conditions have been determined such that the energy per 
baryon is given by energy and baryon flow, and the total baryon number is 
i](Ai + A%), as stopped in the interaction region. They are shown in table |[ 
The radius are for zero impact parameter. For this, equations of state of the 
QGP are needed, and we have employed our modelES in which the perturbative 
correction to the number of degrees of freedom were incorporated along with 
thermal particle masses. n 

Allowing for dilution of the phase space density in expansion, we integrates 
a population equations describing the change in %(t): 
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Table 2: The initial conditions for S-Pb/W at 200A GeV and Pb-Pb at 158A GeV for 
different stopping values r). 





t in [fm] 


V 


Rin [H 


Tl n [MeV] 


A q 


S-Pb/W 


1 


0.35 


3.3 


280 


1.5 




1 


0.50 


3.7 


280 


1.5 


Pb-Pb 


1 


0.50 


4.5 


320 


1.6 




1 


0.75 


5.2 


320 


1.6 



Table 3: 7 S and N s /B in S-W at 200A GeV and Pb-Pb at 158A GeV for different stop- 
ping values of baryonic number and energy t]b = VE computed for strange quark mass 
m s (lGeV) = 200 ± 40 MeV, a s (M z ) = 0.118 . 



-Blab 


S-W at 200A GeV 


Pb-Pb at 158 A GeV 


vb = m 


0.35 0.5 


0.5 0.75 


Is 


0.53 ±0.14 0.65 ±0.15 


0.69 ± 0.15 0.76 ±0.16 


N s /B 


0.67 ± 0.16 0.70 ±0.16 


0.954 ±0.20 0.950 ±0.20 



Here K2 is a Bessel function and x = m s /T. Note that even when 1 — 7^ < 1 
we still can have a positive derivative of 7 S , since the first term on the right 
hand side of Eq. (^7|) is always positive, both T and d/dx(x 2 K2) being always 
negative. This shows that dilution due to expansion effects in principle can 
make the value of 7 S rise above unity. 

Given the relaxation constant T s (T(i)), these equations can be integrated 
numerically, and we can obtain for the two currently explored experimental 
systems the values of the two observables, 7 S and N s /B, which are given in 
table [?]. There is a considerable uncertainty due to the unknown mass of strange 
quarks. However, since this is not a statistical but systematic uncertainty 
depending nearly alone on the value of the strange quark mass parameter, all 
the results presented will shift together. We note further that there seems to be 
very little dependence on the stopping fractions in the yield of strange quarks 
per baryon N s /B. Thus if the expected increase in stopping is confirmed, we 
should also expect a small increase by 15% in specific strangeness yield. 

5 Hadronization from Fireball 
Fireball parameters 

We next introduce all the model parameters used in the fit of the particle ratios, 
not all will be required in different discussions of the experimental data. For 
more details about the thermo-chemical parameters we refer to the extensive 
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discussion in the earlier study of S-S and S-W data 0. The key parameters 
are: 

1) Tfi the formation/emission/freeze-out temperature, depending on the 
reaction model. Tf enters in the fit of abundance ratios of unlike particles 
presented within a fixed p± interval. The temperature Tf can in first approxi- 
mation be related to the observed high-mj_ slope Tj_ by: 

T± ~ T f - 1 + V± -. . (48) 
ll - v\ - vjj 

In the central rapidity region the longitudinal flow vu ~ 0, in order to assure 
symmetry between projectile and target. Thus as long as Tf < Tj_, we shall 
use Eq. ( [f8| ) setting U|| = to estimate the transverse flow velocity v± of the 
source. 

2) A q : the light quark fugacity. We initially used in our fits both u, d- 
flavor fugacities A u and Ad , but we saw that the results were equally adequate 
without allowing for up-down quark asymmetry, using the geometric average 
A q = vAu^d j moreover the fitted up-down quark, fugacity asymmetry was 
found as expected in our earlier analytical studies E3. 

3) A s : the strange quark fugacity. A source in which the carriers of s 
and s quarks are symmetric this parameter should have the value A s ~ 1, in 
general in a re-equilibrated hadronic matter the value of A s can be determined 
requiring strangeness conservation. 

4) 7 S : the strange phase space occupancy. Due to rapid evolution of dense 
hadronic matter it is in general highly unlikely that the total abundance of 
strangeness can follow the rapid change in the conditions of the source, and 
thus in general the phase space will not be showing an overall abundance equi- 
librium corresponding to the momentary conditions. 

5) Rq\ in table ^ we note this parameter describing the relative off- 
equilibrium abundance of strange mesons and baryons, using thermal equi- 
librium abundance as reference. 

Parameter R s c is needed, when we have constraint on the strangeness abun- 
dance and/or when we address the abundance of mesons since there is no 
a priori assurance that the relative production/emission strength of strange 
mesons and baryons should proceed according to relative strength expected 
from thermal equilibrium. Moreover, it is obvious that even if re-equilibration 
of particles in hadronic gas should occur, this parameter will not easily find 
its chemical equilibrium value i?Q = I as we alluded to in section 1. How- 
ever, due to reactions connecting strange with non-strange particles we expect 
Rq = Rq, where Rc is the same ratio for non-strange mesons and baryons, 
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using thermal abundance as reference. The value of Rq > 1 implies meson ex- 
cess abundance per baryon, and thus excess-specific entropy production, also 
expected in presence of color deconfinement L3. 

The relative number of particles of same type emitted at a given instance 
by a hot source is obtained by noting that the probability to find all the j- 
components contained within the i-th emitted particle is 

■v, x-;ll-v ' ' • (49) 

jei 

and we note that the total energy and fugacity of the particle is: 

jei jei 

The strangeness occupancy 7 S enters Eq. ( ff9| ) with power k, which equals 
the number of strange and anti-strange quarks in the hadron i. With Ei — 
\J mf + p 2 = \J rrij + p\ cosh y we integrate over the transverse momentum 
range as given by the experiment (here p± > 0.6 GeV) taking central rapid- 
ity region y ~ to obtain the relative strengths of particles produced. We 
then allow all hadronic resonances to disintegrate in order to obtain the final 
relative multiplicity of 'stable' particles required to form the observed parti- 
cle ratios. This approach allows to compute the relative strengths of strange 
(anti)baryons both in case of surface emission and equilibrium disintegration 
of a particle gas since the phase space occupancies are in both cases properly 



accounted for by Eq. (49). The transverse flow phenomena enter in a similar 
fashion into particles of comparable mass and are not expected to influence 
particle ratios. Finally we note that particles which are easily influenced by 
the medium, such as (j), require a greater effort than this simple model, and 
are also not explored in depth here. 

We obtain the least square fit for the eight above reported (anti)baryon 
ratios. Our first approach is motivated by the reaction picture consisting of 
direct emission from the QGP deconfined fireball. The value of statistical pa- 
rameters controlling the abundances are thus free of constraints arising in an 
equilibrated hadronic gas (HG) state l3. The fitted thermal parameters are 
presented in the first line of table || along with the total x 2 f° r the eight ratios. 
The fit is quite good, the error shown corresponds to the total accumulated 
error from 8 measurements; even if one argues that it involves 4 parameters 
to describe 5 truly independent quantities, the statistical significance is con- 
siderable, considering that 8 different measurements are included. Such a free 
fit does not know that we are expecting that the final state comprises a bal- 
ance (s — s) — 0. In order to estimate what would be implied by strangeness 
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Table 4: Values of fitted statistical parameters within thermal model, for 158A GeV Pb— Pb 
strange particle production data. Superscript star a fixed input value for equilibrium 
hadronic gas; subscript '|c': value is result of the imposed strangeness conservation con- 
straint. x 2 i s the total relative square error of the fit for all data points used. First result 
line: direct emission QGP model, no meson to baryon ratio fit. Second line: same, but with 
strangeness conservation yieldin g A =. and variable. Line three: as in line two, in addition 
the meson to baryon ratio Eq. (B3) is fitted. Line four: hadronic gas fit including the ratio 
Eq. Q. 



T f [MeV] 


K 


A s 


7s 


R a c 


x 2 


272 ± 74 


1.50 ± 0.07 


1.14 ± 0.04 


0.63 ± 0.10 




1.0 


272 ± 74 


1.50 ± 0.08 


1-14 c 


0.63 ± 0.10 


4.21± 1.88 


1.0 


151 ± 10 


1.54 ± 0.08 


1.13| C 


0.91±0.09 


0.85±0.22 


1.5 


155 ± 7 


1.56 ± 0.09 


1-14 c 


1* 


1* 


7.6 



conservation constraint among emitted hadrons we present in second line of 
table H the result of a fit assuming that the value of A s is result of the con- 
servation constraint (s — s) =0, and allowing ^ 1, for there should be no 
chemical equilibrium among the emitted strange mesons and strange baryons 
in a sudden and early QGP disintegration. The statistical error is found the 
same as in line 1, since this approach substitutes one parameter (A s ) by an- 
other {Rq). The implication of this fit is that there must be either an excess 
of strange mesons or depletion of strange baryons compared to thermal equi- 
librium expectation, since as we recall R s c — C'^/Cg with d being the yield 
of particles, normalized to one for thermal equilibrium yield. In any case we 
see that strange meson excess is required, which is consistent with excess of 
entropy production. 

The errors seen in the two first lines of table ^ on the statistical parameters 
arise in part from strong correlations among them. In particular the very large 
error in Tf arise from strong (anti) correlation with 7 S : we find that a change 
in one of these parameters is associated with a change in the other at 80% 
level. Further information about the relation of Tf and 7 S may be garnered 
from theoretical considerations. We evaluate using our dynamical strangeness 
production model in QGP how the value of 7 S depends on the temperature of 
particle production Tf. The most important parameter in such a theoretical 
evaluation is the initial temperature at which the deconfined phaseJs created. 
As noted above, we estimated this temperature at T; n = 320 MeVtEl. Further 
uncertainty of the calculation arises from the strange quark mass taken here 
to be m s (lGeV) = 200 MeV. We recall that the strength of the production 
rate is now sufficiently constrained by the measurement of ot s (Mz)- We choose 
a geometric size which comprises a baryon number B = 300 at A q ~ 1.5, and 
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T. = 320 MeV 
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Figure 5: QGP strangessnes occupancy 7 S as function of temperature Tf at time of particle 
production, for intial temperature T n = 320 MeV, with 7s (T n ) = 0.1. 



have verified that our result will be little dependent on small variations in B. 
We show in Fig.|| how the computed 7 S depends on formation temperature 
Tf. The cross to the right shows our fitted value from line 1 or 2 in table || 
It is gratifying to see that it is consistent with the theoretical expectation for 
early formation of the strange (anti)baryons. The relative smallness of 7 S , 
despite the hight strangeness yield, is clearly related to the high temperature 
of particle production. We will below discuss the opposite scenario: the late 
freeze-out of strange particles. 

The relatively high value of temperature Tf we obtained in the reaction 
picture with primordial particle emission is consistent with the experimental 
inverse slopes of the spectra. We are thus led to the conclusion that as far as 
the fitted temperatures and slopes are concerned, it is possible that the high 
m± strange (anti)baryons we have described could have been emitted directly 
from a primordial (deconfined) phase before it evolves into final state hadrons. 

We now consider at the experimental inverse m± slopes. In the common p± 
range of WA97 and NA49 experiments the transverse mass spectrum of A and A 
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obtained by NA49 is very well describe by an exponentialll3. A thermal model 
motivated fit of the inverse slope (temperature) yields = 284 ± 15 MeV 
and Tj^ = 282 ± 20 MeV. This is consistent with the mid-rapidity proton and 
antiproton slope of the NA44 experiment: T p = 289 ± 7 MeV and T ¥ = 278 ± 9 
MeV. For S + jsLa consistent value T H = 290 MeV is also quoted by the NA49- 
collaboration l3. We note that because the baryon masses are large, all these 
slopes are at relatively high m± > 1.3 GeV (for nucleons, in NA44, m± > 1 
GeV). Systematically smaller inverse-transverse slopes are reported at smaller 
m x , for kaons Tf ~ 213-224 MeV for 0.7 < m x < 1.6 GeV in NA490 and 
Tf + = 234±6, Tf ~ = 235±7 MeV in NA440; and 155-185 MeV for 7r,H0, 
depending on the range of p±, but here we have to remember that pions are 
known to arise primarily from resonance decays. An increase of T with m± is 
most naturally associated with the effects of transverse flow of the source. 

Is QGP primordial emission hypothesis also consistent with the chemical 
fugacities we have obtained? The chemical condition is fixed to about 5% pre- 
cision, and there is 40% anti-correlation between the two fugacities A q and A s . 
The information that A s 7^ 1 is contained in at least two particle abundances; 
arbitrary manipulation of the reported yields of one particle abundance did 
not reduce the value A s to unity. Since A s 7^ 1 by 4s.d. it is highly unlikely 
that A s = 1 is found after more data is studied. While one naively expects 
^QGP = i ( t assure the strangeness balance (s — s) — 0, there must be a 
small deviations from this value, even if the emitted particles were to reach 
asymptotic distances without any further interactions: in presence of baryon 
density the deconfined state is not fully symmetric under interchange of par- 
ticles with antiparticles. A possible mechanism to distinguish the strange and 
anti-strange quarks arises akip_.to the effect considered for the K~ /K + asym- 
metry in baryonic matter E3c3: there is asymmetric scattering strength on 
baryon density which causes presence of a mean effective vector potential 
W . Similarly, strange quark interaction with baryon density would lead to a 
dispersion relation 

E a/S = yfi%+?±W, (51) 

and this requires in the statistical approach that the Fermi distribution for 
strange and anti-strange quarks acquires a compensating fugacity A s § = e^- w l T 
to assure strangeness balance in the deconfined phase. In linear response ap- 
proach W oc Vh consistent with both W and baryon density = (n q — riq)/3, 
being fourth component of a Lorentz-vector. It is clear for intuitive reasons, 
as well as given experimental observations, that the baryon stopping and thus 
density increases considerably comparing the S and Pb induced reactions in 
the energy domain here considered. We also recall that in S-W reactions 



3G 



Ag ~ 1.03 ± 0.05E3. Should in the dense matter fireball the baryon density 
grow by factor 2-4 as the projectile changes from S to Pb, this alone would 
consistently explain the appearance of the value A s = 1.14 ± 0.04 obtain using 
W oc v\, scaling. It is worth noting that the value W ~ 38 MeV suffices here. 
Note also that the Coulomb potential effect on the charge of the strange quarks 
is of opposite magnitude and about 1/5-1/6 of the here required strength. 



Late emission scenario: HG with or without QGP? 

A generally favored picture of particle production involves flow expansion of 
the primordial phase till a transition temperature is reached, at which time 
the final state hadrons are produced, and soon thereafter freeze out. These 
particles may directly reach a detector or re-equilibrate and appear to the 
observer as if emitted from a HG phase, except that entropy/strangeness excess 
effect should remain. In order to force our data fit to converge to such a late 
particle production scenario we introduce the experimental result, which was 
essential for such an argument in the S-Pb induced reactions. The quantity 
of interest is ratio of particle yields between particles of very different mass. 
Thus in lines 3, 4 in table [| we include in the fit also as experimental input 
an estimate of the hyperon to kaon ratio, thus altogether we now fit 9 data 
points. _ 

We note that the NA49 spectra ej of kaons and hyperons have a slightly 
overlapping domain of m±. We recall that the slopes are not exactly equal, 
thus all we can do is to try to combine the two shapes, assuming continuity 
consistent with flow, and to estimate the relative normalization of both that 
would place all experimental points on a common curve. We have carried out 
this procedure and obtained: 



_A_ 



6.2 ±1.5. (52) 



Note that there is a tacit presumption in our approach that a similar effective 
Ay interval wasjised in both spectra. We recall that this ratio was 4.5 ± 0.2 
in the S-W dataS. 

Our approach in the third line corresponds to a freeze-out from 'cold' QGP 
phase, in that we allow the abundance parameters 7 S and Rq to deviate from 
HG equilibrium values. We note that this cold-QGP alternative has a very 
comparable statistical significance as the hot-QGP. Given the low tempera- 
ture and high m± inverse slopes we must have considerable transverse flow. 
The computed flow velocity at freeze-out is Vf = 0.51 . This is just below the 
relativistic sound velocity v s = 1/V3 = 0.58, which we have assumed. In 
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Fig. | the cross to the left shows the result of the fit we just described; al- 
lowing for potentially smaller expansion velocity and all the above discussed 
uncertainties in the computation, this result must also be seen as a very good 
agreement between the result of data fit and the theoretical calculation. This 
also means that we cannot distinguish in the present data between early forma- 
tion of strange antibaryons and an expansion model followed by direct global 
hadronization. 

Since the fitted values of 7 S and Rq allow the HG equilibrium, we attempt 
such a fit in line 4, where the particle yields are fitted constrained for HG 
equilibrium, and we use strangeness conservation to evaluate the strangeness 
fugacity A s . We show the result of the fit in the last line of table [| and in 
particular we note: 

T f = 155 ±7 MeV, -> v± ~ 0.5 ~ v s ; 
A q = 1.56 ±0.09, -► A s = 1.14; 
X 2 /9 =0.84, -► C.L. >60%. (53) 

We recall that the baryo-chemical potential is given in terms of T and A q , 
specifically /Xb = 3Tln A q , and we find /ib = 204 ± 10 MeV in this hadronic gas 
condition. We note that while the quality of the fit has degraded, it still has 
considerable statistical significance. It would appear to be a 'good' fit to the 
naked, an unequipped eye. 

Unlike S-induced reactions, the hypothesis of an equilibrated HG in the 
final state cannot be easily discarded in Pb-Pb collisions since the thermo- 
chemical parameters are extraordinarily consistent with this hypothesis and 
the principle of strangeness conservation, as we illustrate in Fig.^. Here, the 
cross corresponds to the fitted properties of the particle source, while the lines 
correspond to the constraint of the HG gas source to yield (s — s) = at 
finite baryon density represented by the value of A q . Thus the strange quark 
fugacity is in general not zero and the cross falls just on freeze-out at T = 160 
MeV when the meson-baryon abundance is equilibrated. The slight difference 
in freeze-out value to results shown in last line of table [| is result of using the 
WA97 A/ A value Eq. @ rather than the averred value Eq. @. 

If HG is indeed present in the final state, the proper interpretation of these 
data, and the likely_reaction scenario, compatible with our earlier work on S 
induced reactions EH is as follows; the relatively large fireball of dense and 
deconfined matter disintegrates and produces dense, confined hadronic gas in 
which strange particles have time to re-scatter and to establish relative chem- 
ical equilibrium. A possible test of this hypothesis would be to see variation 
of the chemical parameters as the size of the fireball changes with impact pa- 
rameter (centrality of collision) since re-equilibration should diminish for small 
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Figure 6: Strangeness conservation constraint in HG as function of freeze-out A q : the lines 
correspond to different freeze-out temperatures Tf (solid 160 MeV, dashed 140 MeV and 
dotted 120 MeV). The cross corresponds to the chemical freeze-out we determined above. 

reaction volume. However, such data are presently not available, and there is 
no indication that indeed a change of the strange (anti)baryon yields occurs 
as the centrality of the interaction is reduced. On the other hand, the specific 
entropy and strangeness should comprise a signal of some new physics should 
formation and expansion of QGP phase, followed by re-equilibration into HG 
phase, and freeze-out have occurred. We will now consider the magnitude of 
these effects: 

Strangeness re-equilibration 

When HG emerges from initial dense QGP phase, the number of strange quark 
pairs does not change, but the phase space density of strangeness changes, since 
the phases are different. Because the HG phase has generally a smaller phase 
space density of strangeness than QGP, to conserve strangeness, there will be a 
jump in the phase space occupancy 7 S during the transformation of QGP into 
HG, as there is a jump in the strange quark fugacity. The important point is 
that this could lead to significantly overpopulated HG phase (7 S > 1). This 
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phenomenon can be easily quantified as follows: the observed value of 7g is 
related to the pre-phase change value 7° by introducing the enhancement 
factor we wish to determine: 

7 P^F 77 Q GP (54) 

A simple way to compute the value of the saturation enhancement factor F 1 
is to study the abundance of strangeness per baryon number before and after 
phase transition. 



a/b\ QGP = 7 QGP 
s/b\ H G 7 HG 



F~f = 1L jjp L - = -^W/m , A q , 7s HG ) ■ (55) 



The last expression arises as follows: on the QGP side the abundance of 
strangeness is to a good approximation proportional to 7°- GP and is the integral 
of the strange quark phase space, we evaluate it assuming that to s /Iqgp — 1< 
There is no dependence on chemical properties of the plasma. On HG side, 
at freeze-out we have to evaluate the strangeness-abundance from the strange 
particle partition function given in Eq. (16) of Ea, supplemented by the now 
relevant term comprising ss-rj, rf, <f> states and their resonances. The sum in- 
cludes a terms proportional to ( 7 ^ G )™, with n=l, 2, 3, indicating strangeness 
content of hadrons. The leading kaon and hyperon term is proportional to 7 ^ G 



and hence we have above result, Eq. (p5|). We thus obtain, combining Eqs. (54) 
andEq. @, 

7s b 



F 2 = F _ j_ 

7 7 7s QGP 7s& 



QGP 



(56) 

HG 



where the right hand side now compares the properties of the two phases at the 
boundary between them and we can evaluate it using the theoretical equations 
of state. In analyzes of an experiment we would take the freeze-out parameters 
determined by the fit to data. 

We show, in Fig.^, the strangeness enhancement factor as function of A q 
for several freeze-out temperatures If = 160, 140, 120 MeV, with A s fixed by 
strangeness conservation constraint. We see that F 1 varies typically between 
1.5 and 2, and is specifically 1.6 for the parameter rage of Pb-Pb collisions 
here discussed. This means that observing the value 7 S ~ 1 really means an 
underlying value 7°- GP — 0.6. Conversely, should we be able to create a longer 
lived or hotter QGP state we could expect to observe in the HG phase 7^ 
as large as 1.5-2. Such over-saturation of the phase space would be a rather 
strong smoking gun pointing to the formation of the QGP phase. 
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Figure 7: Strangeness phase space enhancement factor as function of HG freeze out A q 
for Tf = 160 (solid line), T t = 140 (dashed line), T f = 120 (dotted line). Computed for 
7^ G = 1, "Is/Tqgp = 1> including in HG phase kaons, hyperons, cascades, r), <p, f2, and 
imposing strangeness conservation constraint to determine A s . 

Entropy and particle excess 

Another way to argue for the formation of QGP in early stages of an expansion 
scenario of the fireball is to measure the specific entropy experimentally, for 
example by measuring the quantity 

(dN+ dN-\ lfdN+ dN-\ , . 

which we have shown to be a good measure of the entropy content 0. We 
note that in the numerator of Dq the charge of particle pairs produced cancels 
and hence this value is effectively a measure of the baryon number, but there 
is a significant correction arising from the presence of strange particles. The 
denominator is a measure of the total multiplicity — its value is different before 
or after disintegration of the produced unstable hadronic resonances. Using as 
input the distribution of— final state particles as generated within the hadron 
gas final state it is founded that Dq ■ S/ B is nearly independent of the thermal 
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Figure 8: QGP (thick line) and HG (thin line, T = 155 MeV) entropy per baryon S/B 
as function of light quark fugacity A q . Dotted lines guide the eye for the here interesting 
values. 

parameters and varies between 4.8, before disintegration of the resonances, to 
3 after disintegration. 

To obtain a measure of the particle excess, we show in Fig.|8| the specific 
entropy per baryon S/B content in dense hadronic matter as a function of light 
quark fugacity A q . The thick line addresses the deconfined QGP phase, the thin 
line the confined HG phase at T — 155 MeV, with the strange quark fugacity A s 
being determined from the strangeness conservation condition. While the QGP 
result is largely independent of temperature (only other, aside of T dimensioned 
quantity, is m s ), the HG result involves the values of all hadron masses and 
hence is dependent on T . The 12 units of entropy difference between the two 
phases for the here interesting range of fugacity A q = 1.5-1.6 implies that we 
should expect an excess of about 3 mesons per baryon if the deconfined phase 
is formed. We compare with a HG at T — 155 MeV; should the HG phase 
of interest be hotter, this difference between QGP and HG grows, since the 
baryon density in HG grows much faster than entropy for the baryon mass 
is well below the temperature range and thus a change in the factor ton/T 
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matters, while the change in m^/T is immaterial. In other work HG phase at 
190 MeV is often considered, and there the difference between QGP and HG 
properties turns out to be as large as factor twoE3. 

6 Discussion and Outlook: QGP 
Vacuum and the QCD energy scale 

We have relegated to the vacuum properties all scales related to hadronic struc- 
ture. But where-from are these scales derived? It seems that at present there is 
indeed no 'low energy' practical understanding of the hadronic scale. While or- 
dinary QED matter scale is governed by the atomic Bohr radius do = a ( T 1 /m e , 
there is in QCD no obvious relation between quark masses and hadronic spec- 
tra and sizes. The scale of the vacuum condensates (see Section ||), and thus 
of the vaOjUim condensation energy which results e.g. from the Lattice-QCD 
approach lIO is a measured property of the physical ground state, and does 
not relate directly to any fundamental property of the QCD action alone; the 
Ao-parameter of the strong interaction coupling constant a s (see Eq. ( |36| ) and 
the middle-section of Fig.|^) has no further meaning, since in fact the strength 
of this coupling at our energy scale derives from an initial value problem: for a 
given measured value of a s at some energy scale, we integrate the renormaliza- 
tion group equations towards the low energy domain (see Section ^) to obtain 
the needed interaction strength, and the definition of Ao-parameter than arises 
in terms of the initial values. 

It is not the Ao-parameter in the strong coupling constant a s but the 
initial point, thus e.g. the grand unification scale, which is defining for us the 
hadronic energy scale. It is quite possible, though not obvious at all, that 
Ao establishes the strength of vacuum fluctuations, and thus one would be 
tempted to conclude that hadronic mass scales are the result of a 'dimensional 
transmutation' mechanism which transports and reduces the grand unification 
scale to our physical domain. This conclusion is not dependent on the scheme 
of the unification model employed, so the reader who is not happy with initial 
value a s = a e and the associated huge scales 10 16 -10 18 GeV can equally well 
use as the input-scale to strong interactions any other initial point, as long 
as the result yields the now well established experimental value c£l a s (Mz — 
0.118 ±0.03. It is important here to stress that there is no scenario that would 
yield the light quark mass as a significant scale, since the scale properties of 
the physical vacuum always dominate. There is moreover no obvious means to 
use the masses m s , m c which would be easy to associate with QCD condensate 
scales, since in such an approach there would be no understanding why just 
two 'medium' and not the other quarks determine the scales, nor is there a 
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noticeable regularity among the mass parameters of matter (quarks, leptons) 
at any scale Ej. 



Nuclear relativistic collisions and strangeness 

Considering these strong and unexpected conclusions about the fundamen- 
tal nature of strong (nuclear) interactions and its far reaching implications, 
we must assure that the foundations on which these arguments are build are 
exceptionally solid. The unique tool to study the QCD vacuum is generally 
accepted to be the process of its change, the 'micro' bang in the laboratory. In 
the relativistic heavy ion collision experiments we seek to prove by finding a 
way of freeing the color charge of quarks that the effect of confinement arises 
from quantum fluctuations in the QCD vacuum structure. We hope and expect 
to produce in these collisions regions of space filled with mobile color charges 
of quarks and gluons, forming 'QGP', the 5th phase of matter comprising a 
Quark-Gluon (color charge deconfined) Plasma. 

Strange particle production is today appreciated as one of the most in- 
teresting hadronic observables of dense, strongly interacting matter and much 
of the current theoretical and experimental effort in study of relativistic nu- 
clear collisions is devoted to this topic. Our work concentrates on the explo- 
ration of the pattern of production and evolution of hadronic particles carrying 
strangeness flavor. Our continued interest in the subject arises from the re- 
alization that the experimentally observed anomalous production of (strange) 
antibaryons cannot be interpreted without introduction of some new physi- 
cal phenomena, and the hope and expectation is that when more systematic 
experimental data is available, we will be able to make the argument for de- 
confincmcnt which will be generally accepted. 

In our work we assume that thermal quark-gluon degrees of freedom are 
at the basis of many of the hadronic particle production phenomena in rela- 
tivistic hadron reactions. Many simple, but subtle experimental observations 
point in this natural direction. For example, many of the measured m± spectra 
in S- and Pb-induced reactions have the same shape for strange baryons and 
antibaryons of the same kind, and even different kinds of particles when con- 
sidered in same range of m± show the same inverse transverse slope (=spectral 
temperature). 

In passing we note that there has been much discussion of the fact that 
depending on particle mass, the slopes of the particle spectra show different 
temperatures. This behavior is possibly arising from the fact that for different 
particles different ranges of mj_ were considered. This is thus evidence for 
significant expansion flow phenomena governing final stages of QGP disinte- 
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gration and in no way does this behavior of the spectra contradict the picture 
of local thermal equilibrium. 

We have recently improved the two particle collision cross sections leading 
to strangeness and charm production. This was done using the recent advances 
in measurement of the strong coupling constant a s at the scale of Mz, a 
result which eliminated the need for the introduction of an arbitrary coupling 
constant in flavor production calculations. Using simple dynamical models 
of dense matter evolution we are able to obtain the final state abundance 
and phase space occupancy of strangeness originating in the deconfined phase. 
From this condition at the point that deconfined phase ceases to exist and 
hadronizes, we develop using conservation laws, constrains on the hadronic 
particle abundances. 

Relative abundances of strange hadrons allowed us to investigate the chem- 
ical equilibrium. We find that there is enough time in the collision which passes 
through the deconfincment phase to establish near-chemical equilibrium con- 
ditions locally. The near chemical equilibrium leads to abundance anomalies 
such as the predicted E3 A/p > 1, found in all A-A collisions eJ. There is 
a priori no opportunity in terms of reaction mechanisms and time constants 
for a confined HG state to reach conditions consistent with this result. This 
result is generally viewed to he-.a_signature for the primordial QGP phase, as 
such abundances are expectedE3l!3 if the HG state of matter has as source the 
deconfined soup of quarks and gluons. 

We did not discuss in depth the entropy excess (see, however, our earlier 
report!), which also are inconsistent with pure HG state, but are well under- 
stood in terms of deconfinement. 

Our view is that the observed hadronic particle effects, and specifically 
strangeness, provide a simple, consistent interpretation of the heavy ion col- 
lision data within the hypothesis that a novel type of deconfined hadronic 
matter is formed in both S- and Pb-induced reactions on heavy nuclei, but 
that quite different initial conditions are reached in these two cases. Moreover, 
the longitudinal flow, clearly visible in the S-S 200A GeV, in e.g., A-rapidity 
spectra, is not as pronounced in Pb-Pb reactions, but there seems to be rela- 
tively strong transverse/radial flow driven by the high internal pressure which 
has been reached in the initially highly compressed fireball matter. We also ex- 
pect that the much greater volume of the fireball formed in Pb-Pb reactions, 
as compared to S-Pb leads to a greater volume of deconfined initial phase, 
which in turn in the hadronization process should enhance the re-equilibration 
of final state hadrons after freeze-out. Thus the Pb-Pb collisions should more 
approach the HG type chemical equilibrium of particle yields than S-W/Pb or 
S-S reactions. 
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The only potentially contradiction to this unified view of the 160-200 GeV 
A heavy ion reactions may be the systematics of charmonium production which 
suggest a difference in the absorption rate of cc in central dense matter formed 
in Pb-Pb interactions as compared to the earlier studies of pA and SA colli- 
sions. There is no full and consistent theoretical understanding of all experi- 
mental details which may yet turn out to confirm our views about the presence 
of deconfinement in all these conditions. 

Immediate future 

Underlying the strange particle signatures of deconfinement is the rather rapid 
chemical equilibration of quark flavor, originating in gluon fusion reactions. 
Gluons themselves are the QGP fraction which is generally believed to most 
rapidly approach chemical equilibrium in the deconfined phase. We presumed 
in our studies that at SPS energies that gluon chemical equilibrium has been 
reached, considering the lifespan of the reaction. Such an assumption is less 
compelling at RHIC energies, indeed some model calculations exist arguing 
that even gluons could not reach full equilibrium. These conclusions are 
reached assuming a fixed and small QCD interaction strength a s . Since there 
is no threshold to production of gluons, and Bose effects enhance in dense mat- 
ter gluon formation, our study of the strong scale dependence of a s suggest 
that one has to review the gluon chemical equilibration in the near future. In 
anticipation of this, we recommend use of glue chemical equilibrium for RHIC 
energy scale collisions as well. 

The dynamics of evolving and exploding QGP phase impacts significantly 
the resulting entropy distribution in rapidity and strange particle abundances. 
At present we have also engaged in several approaches to the explosive flow 
problem: we study 3-d hydrodynamic solutions, develop semi-analytically sol- 
uble homologous hydrodynamic evolution/flow models of QGP pad study the 
entropy sources during the flow phase. Our schematic approaches will permit 
by its simplicity a systematic investigation of the complex dynamics of strongly 
interacting hadronic matter evolution comprises several interwoven challenges: 

a) the determination of initial flow conditions; 

b) description of explosive flow of dense matter and the radiation reaction; 

c) influence on matter flow by the equations of state 

d) constrains on the final state conserved quantities (baryon number, entropy, 
and strangeness). 

While hydrodynamical flow of thermally equilibrated matter is not entropy pro- 
ducing, other dynamical processes occur during the evolution of dense matter 
that can produce entropy, for example chemical equilibration, i.e. approach 
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of particle abundances to their equilibrium number is entropy producing. We 
consider of considerable importance that relativistic hydrodynamical dynami- 
cal flow equations are improved to allow for entropy production by approach 
to phase space saturation by gluons and quarks. 

The key problem we need to resolve is how to establish that we have indeed 
seen the thermal, deconfined phase, where the degrees of freedom are those of 
nearly massless quarks and gluons? Other forms and phases of hadronic matter 
can perhaps yet be invented that could produce similar observational effects. 
We need to devise experimental approaches that would allow us to understand 
the rich deconfinement properties of dense relativistic nuclear/quark matter. 
In immediate future we have to explore the collision energy dependence and see 
how the key properties change between the low AGS energies, the intermediate 
SPS energies towards the immediate future of RHIC. 
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